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Abstract. We analyse the nucleation of quark matter droplets under protoneutron 
star conditions. We adopt a two-phase framework in which the hadronic phase is 
described through a non-linear Walecka model and the just deconfined matter by 
the MIT bag model including color superconductivity. Surface tension and curvature 
energy are calculated self-consistently within the MRE formalism. We impose flavour 
conservation during the transition, which means that the just deconfined quark droplet 
is transiently out of equilibrium with respect to weak interactions. Our results show 
that trapped neutrinos slightly increase the critical density for deconfinement and that 
color superconductivity significantly decreases such density at low temperatures. We 
also show that the nucleation rate is negligible for droplets larger than 100 — 200 fm 
and is huge for smaller droplets provided that the temperature is low enough. We 
compare our results with previous calculations using the Nambu-Jona-Lasinio model 
with color superconductivity and the MIT bag model without color superconductivity. 
We conclude that the deconfinement transition should be triggered instantaneously 
when a density slightly larger than the bulk transition density is reached at some layer 
of a protoneutron star. Since color superconductivity lowers the transition density at 
low temperatures, the transition is likely to occur after significant cooling in a massive 
enough protoneutron star. 
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1. Introduction 

Understanding the behaviour of the strong interacting matter for moderate temperatures 
and high baryonic density is one of the main tasks in astrophysics for describing the 
interior of neutron stars and regions close to the core of collapsing stars. In the interior 
of that objects, at high densities, the deconfinement transition from hadronic to quark 
matter might occur. This transition, in astrophysical compact and dense objects, is 
basically studied in two scenarios: in protoneutron stars (PNS) and neutron stars (NS). 
The PNS are compact objects remaining after the gravitational collapse and supernovae 
explosion of a massive star. Initially the temperature is about few tens of MeV and the 
neutrinos are temporarily trapped in the interior of the PNS. After a few seconds of 
evolution, the neutrinos are radiated, and the PNS cools down turning into a NS. The 
interior of these objects can reach a density well above the nuclear saturation density 
(po), and in this scenario the hadronic matter can suffer a deconfinement transition into 
quark matter [U El [31 HI [5] . The phase transition should begin with the nucleation of a 
small (~ 10~ 14 m) quark-matter droplet near the center of the star [6], d El [9], [TUl [TT] 
that may later grow through a combustion process converting a macroscopic portion of 
the star into quark matter [12j . 

Even though hadrons are a bound state of quarks, we can not describe with a single 
model the thermodynamics of the transition from hadronic to quark matter. Then, we 
need to adopt different models to describe each phase. In the present paper we analyse 
the deconfinement transition in protoneutron star conditions employing the MIT bag 
model in the description of quark matter. For the hadronic phase we use a model based 
on a relativistic Lagrangian of hadrons interacting via the exchange of a, p, and u 
mesons. The nucleation is treated as a first order phase transition and the finite size 
effects of creating a droplet of deconfined quark matter in the hadronic environment are 
described using the multiple reflection expansion (MRE) framework [131 El El El E] ■ 
Two important features we consider in our work are: (a) quark flavor is conserved during 
the deconfinement transition because it is driven by strong interactions [6l El El |9l El E] 
and (b) when color superconductivity is included together with flavor conservation, the 
most likely configuration of the just deconfined phase is two-flavor color superconductor 
(2SC) provided the pairing gap is large enough pO, [13]. 

The article is organized as follows. In Sec. 2 we present the main aspects of the 
non-linear Walecka model describing the hadronic phase. In Sec. 3 we present the 
generalities of the model we use for the quark phase and the MRE formalism. In Sec. 
4 we study the deconfinement transition at finite temperature for different neutrino 
trapping conditions. In Sec. 5 we present our results followed by a summary and 
conclusions in Sec. 6. 

| Works that deal with the structure of hybrid proto- neutron stars [18) . consider the hadron-quark 
interphase as being in equilibrium under weak interactions. This condition is appropriate for such 
studies but it isn't for the present analysis of the deconfinement of quark droplets. 
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2. Hadronic matter 

We use a non-linear Walecka model for describing the hadronic phase (see e.g. [3] and 
references therein) composed by the baryon octet (n, p, A, S + , S°, E~, H~, 5°), electrons 
e~, electron neutrinos v e and the corresponding antiparticles. The Lagrangian for this 
model is the following [19J 



= ^ B ^{idi* ~ XuB g^Uy. ~ x p b 9 p t- py) 

B=n,p,A,S+.°.-,H-. 



- [m B - x aB g a <y) 



1 be 
ipB + -(d^o-d^a - m\o 2 ) - -m N (g a a) 3 - -{g^af 



-u^" + -rrfor - -p> • + - m; p M • p< 
+ Yl Mill*® 1 * -m L \il) L , (1) 

L=e,v e 

where B, M and L refer to baryons, mesons and leptons respectively. Baryons interact 
through the exchange of mesons a, co and p. The constants x a B = gaB/da, x w b = QwbI ' 9uj 
and x p b = g P B/ 'g P are the ratios of the coupling constants of the hyperons to the coupling 
constants of the nucleons. The EoS is obtained from the above Lagrangian by means of a 
relativistic mean field model treatment (see [3] and references therein for more details). 
In order to adapt the equation of state (EoS) to the actual conditions prevailing in 
protoneutron stars, the lepton term includes the contribution of electron neutrinos. 
The thermodynamic potential for this model is: 

2 , / \ 2 



B,L 

+ i6m n (, CT a) 3 + ic(^) 4 -i(^y P ^. (2) 

The weighted isospin density p' l3 and the weighted baryon density p B are given by: 

Pu = /^pihirij, (3) 

i=B 

pB = ^2 x ^n i} (4) 

i=B 

where I^i is the third component of the isospin of each baryon and rij is the particle 
number density of each baryon: 

* = j£yj UiiT) - UT% (5) 

being 7, the degeneracy factor. 

The mean field g a a satisfies the equation: 

2 

( — ) (9a0-) + bm n (g a a) 2 + c{g a af = ^a^nf, (6) 

771(7 J i=B 
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where nf is the scalar density: 

From the thermodynamic potential we obtain the total pressure P and the mass - 
energy density p: 

1 / 9u \ 2 ' 2 1 ( 9a N - 1 



E p i + o(—) t>&-^[—) M 



i=B,L 



P «= Ep , + i(j^)Vi + i(^)-W 

2\m„,J 2\mJ 



i=B,L 



+ \bm n {g c af + \c(g a aY + & (9) 

Here Pj and pi are the expressions for a Fermi gas of relativistic, non-interacting particles: 

p -=\i^l d3p w^w (/ - (r)+/ - (r)) - (10) 

ft = / <?J> (P 2 + m* 2 )" 2 (/.(T) + /,(T)), (11) 

where fi(T) and /j(T) are the Fermi— Dirac distribution functions for particles and 
antiparticles respectively: 

U{T) = (exp([(p 2 + ruff - p*]/T) + (12) 

UT) = (exp([(p 2 + mf ) 1 / 2 + p*]/T) + (13) 

Note that for baryons we use, instead of masses mj and chemical potentials pi, "effective" 
masses m* and chemical potentials p* given by: 

m* = rm- x ci {g c a), (14) 

p* = pi- x ui (—) p B - Xpihi (^-) p'j (15) 
Vm^/ \m p J 6 

The hadronic phase is assumed to be charge neutral and in chemical equilibrium 
under weak interactions. Electric charge neutrality states: 

n p + ns+ — ri£- — ra s - — n e = 0. (16) 

Chemical weak equilibrium in the presence of trapped electron neutrinos implies that 
the chemical potential pi of each baryon in the hadronic phase is given by: 

A*i = QBl^n ~ q e (Ve ~ AO» (17) 

where qs is its baryon charge and q e is its electric charge. For simplicity we are assuming 
that muon and tau neutrinos are not present in the system, and their chemical potentials 
are set to zero. 



Nucleation rate of color superconducting droplets in protoneutron stars 



5 



The values of the five constants of the model are determined by the properties of 
nuclear matter. Three of them determine the nucleon couplings to the scalar, vector 
and vector-isovector mesons g a /m ai g^/m^, g p /m pi and the other two determine the 
scalar self interactions h and c. Moreover, we consider that all hyperons in the octet 
have the same coupling as the A, and that the coupling of the A is 90 % of that of the 
nucleons. Therefore, we have that x a B, x uB and x p b are equal to 1 for the nucleons and 
0.9 for hyperons. In this paper we use the parametrization labelled as GM4 in [13] with 
the following values: (g a /m a ) 2 = 11.79 fm 2 , (g w /m w ) 2 = 7.149 fm 2 , (g p /m p ) 2 = 4.411 
fm 2 , b = 0.002947, c = -0.001070. With this parametrization the EoS is stiff and gives 
a maximum mass of 2.0 M & for compact stars. 

All the above equations can be solved numerically by specifying three 
thermodynamic quantities, e.g. the temperature T, the mass-energy density p H and 
the chemical potential of electron neutrinos in the hadronic phase ^ . 



3. Quark matter 

3.1. Formalism for quark matter in bulk 

The quark phase is composed by u, d, and s quarks, electrons, electron neutrinos and 
the corresponding antiparticles. We describe this phase using the MIT bag model at 
finite temperature with zero strong coupling constant, m u ^ = and strange quark mass 
m s = 150 MeV. The total thermodynamic potential for the quark matter phase can be 
written as 



E 



n 



cf 



(18) 



where / = u, d, s is the flavor index, c = r,g,b is the color index and L stands for the 
leptons. The contribution for free unpaired quarks is given by 

7 T 



n 



cf 



2vr 2 



r In 



1 + e 



dk, 



(19) 



being E c f = \/k 2 + m 2 c f the particle energy, and /i c / the chemical potential. In the case 



of paired quarks we consider [20 

7 T 



n 



2tt 2 



/•oo 

/ k 2 In 


1 + e t 


/o 





dk, 



(20) 



where s c f = ±a/ (E c f — /i c /) 2 + A 2 is the single-particle energy dispersion relation when 
it acquires an energy gap A. Note that, for particles, we can obtain Eq. ([19]) from Eq. 
(120p in the limit A = 0, by taking the minus sign in the dispersion relation for E c f < fi c f 
and the plus sign when E c f > fi c f (see e.g. |20j). 

The gap equations for a color-superconducting condensate with total spin J = 
have been derived perturbatively in dense QCD [21] . To leading order in weak coupling, 
the temperature dependence of the condensate is identical to that in BCS-like theories 
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. Then, we consider the following temperature dependence for the gap parameter in 
Eq. (J20D 




A(T) = A ot /l-(^) , (21) 



where the critical temperature is T c = 0.57A [21 J . 

Each lepton species contributes with a term of the form 



l + e 



dk, (22) 



with L = e~, e + , u e , u e and E L = \Jk 2 + m\. The degeneracy factor is 7 = 2,2, 1 for 
quarks, electrons and neutrinos, respectively. In all cases, the antiparticles contribution 
is obtained considering p, = —fi. 

In the present model, the quantities B, m s and Ao are free parameters. According 
to their values, the energy per baryon of three flavor deconfmed matter (composed 
of quarks u, d and s) at zero pressure and temperature can be higher or lower than 
the mass of the neutron m n . This condition defines the so called stability windows 
[22] [23] . which are the regions in the m s — B parameter space where quark matter 
is self-bound, i.e. it is the true ground state of strongly interacting matter. In this 
paper we use B = 100 MeV /fm 3 corresponding to absolutely stable quark matter and 
-B[MeV/fm 3 ] = 160,353 corresponding to quark matter allowed only at high pressures. 
The case B = 353 MeV/fm 3 leads to an equation of state rather similar to set 1 of the 
N JL model used in Refs. [U [13] . 

The here-considered just deconfined phase is out of chemical equilibrium with 
respect to weak interactions. As we will see in the next section, the chemical potentials 
fi c f, fi e and fi Ue are related among each other through flavour conservation, color 
neutrality and pairing conditions. 

3.2. Finite size effects 

In order to study the formation of finite size droplets of quark matter we use the multiple 
reflection expansion (MRE) formalism [131 El E3 HEl EE] • For a finite spherical droplet, 
the modified density of states is given by 

6n 2 127r 2 

PMRE(k,m f ,R) = 1 + —f s + jj-^fc (23) 



with 



2 / k \ 

1 arctan — ) 

7r \mj / 



(24) 



fc = T^ 1-7^^-arctanfAV) . (25 ) 



8tt 

1 

12~7T 2 r 2m f V 2 \mfJJl 
The fs and fc terms correspond to the contributions of surface and curvature, 
respectively, and the drop radius is given by R. As in [13], we use the Madsen ansatz 
[T5] to take into account the finite quark mass contributions. 
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The MRE formalism is included in the equation of state by means of the following 
replacement in the thermodynamic integrals 



d k ->■ / d 6 k p MRE . (26) 
o Ja ir 

According to Ref. [HUE], the density of states for massive quarks is reduced if compared 
with the bulk density, leading to a negative density of states for a range of small 
momentum. In order to avoid this nonphysical effect one can introduce an infrared 
cutoff A IR in momentum space k |T6l IT7] . For calculating the value of A IR , we need to 
solve the equation pmre = with respect to the momentum k and take the larger root 
as the IR cutoff. Then, this cutoff depends as well on the radius of the spherical drop 
and the quark mass. 

Performing the above replacement into Eq. (|T9l) . the thermodynamic potential for 
the quark matter phase reads 

Q MRE = _pQyQ + aA + ( - C (27) 

Here, A = 4nR 2 is the spherical drop area, C = 8nR the curvature and V® = §7ri?' 
the volume. The pressure P® is given by 

8Qq RE 



3 



pQ 



the surface tension is 



(2* 



a = 



OA 



T,!i,VQ,C 



and the curvature energy density is obtained as follows 



C 



DC 



T,v,A,VQ 



(29) 



(30) 



4. Nucleation of quark droplets in hadronic matter 

To calculate the conditions for the transition we assume thermal, chemical and 
mechanical equilibrium between the quark droplet and its hadronic environment. 
Thermal and chemical equilibrium means that the temperature T and the Gibbs free 
energy per baryon g are equal in both the hadronic (H) and the quark matter (Q) 
phases: 

t h =T Q = T, (31) 

g H (T,p H ,^J =9 Q (T,{^c},^,^J- (32) 
To obtain the condition for mechanical equilibrium we write the total thermodynamic 
potential for the combination of a just nucleated quark matter drop immersed in a 
homogeneous hadronic environment: Q = Qh + Qq RE ', where Qh = —P H V H for the 
hadronic phase. The condition for mechanical equilibrium is obtained from dQ/dR = 



Nucleation rate of color superconducting droplets in protoneutron stars 



8 



and reads 



11 



P Q (T,{ Hc },^,^J-P H (T,p H ,^ 
!<r(r, 0* /c }, a*?, /£) + ^C(T, {/i /c }, a*?, /£)■ (33) 

On the other hand, deconfinement is driven by strong interactions and therefore 
quark and lepton flavours must be conserved during the deconfinement transition 
P d El El EH [TT] . When a small quark- matter drop is nucleated at the core of a 
compact star, the abundances of all particle species inside it must be initially the same 
as in the hadronic matter from which it has been originated. Thus we have 

Y f H (T,p H ,^J = ^ CM/Vc}, (34) 

with / = u,d, s,e,u e , being Yj 1 = nf/n% and Y9 = nj /n% the abundances of 
each particle species in the hadronic and quark phase respectively. Notice that, since 
the hadronic phase is assumed to be electrically neutral, flavor conservation ensures 
automatically the charge neutrality of the just deconfined quark phase. As a consequence 
of Eq. (l34"j) . in the present work, the just deconfined phase is out of (3 equilibrium. 

Additionally, the deconfined phase must be locally colorless, i.e. it must be 
composed by an equal number or red, green and blue quarks: 

n r (T, {n fr }) = n g (T, {n fg }), (35) 
n r (T, {fifr}) = n b (T, {///&}). (36) 

Finally, it has been shown that when color superconductivity is included together 
with flavor conservation and color neutrality, the most likely configuration of the just 
deconfined phase is 2SC provided the pairing gap is large enough [7J. Thus, in order to 
allow for pairing between quarks d r with u g and between quarks u r with d g we impose 
that: 

n ur (T,/i ur ) = n dg (T,jj, dg ), (37) 
n dr (T, n dr ) = n ug (T, fi ug ). (38) 

Equations fl31H38p together with the assumption that p sr = fj, sg = fi s b allow to determine 
the transition density for given values of the droplet radius, the temperature and the 
chemical potential of the neutrinos in the hadronic phase. 

The theory of homogeneous nucleation predicts a rate of nucleation of a droplet of 
radius R as given by [T5l [24] 

r « T 4 e- AQ/T , (39) 

where AQ, the free energy involved in the formation of the droplet, is given by 

ACl = -^ttR 3 (P q - P H ) + 47T-RV + 8nR(. (40) 

The prefactor T 4 is included in Eq. ( 139|) on dimensional grounds because the 
nucleation rate is largely dominated by the exponential [H [T31 1221 ES], i-e. we always 
have log lo r w log 10 (prefactor) — Afi/[Tln(10)] with the second term much larger than 
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the first. A more elaborate statistical prefactor has been developed in the literature 
[27] . However, the determination of such prefactor involves the specification of transport 
coefficients such as the thermal conductivity, and the shear and bulk viscosities, which 
are not well known for ultradense matter (see Refs. [TJ [27] for more details). We have 
checked that the prefactor given in Ref. pQ is in fact very different from the T 4 factor, 
but it is not dominant with respect to the exponential for the conditions encountered 
in our calculations. Thus, our results are not significantly affected by the choice in Eq. 
( 1391 . The nucleation time r is the typical time needed to nucleate a droplet of radius 
R and is given by r = (^ttR 3 

5. Results 

In Fig. [Tjwe show the mass-energy density and the temperature at which the conversion 
of a portion of hadronic matter into a quark droplet is energetically favourable. By 
comparing left and right panels, we can see that curves with /z^ = 200 MeV are displaced 
to the right with respect to the same curves for /i^ = MeV. In other words, the trapped 
neutrinos push the transition density to higher values, i.e. the trapped neutrinos tend 
to inhibit the transition, although the effect is not very large. This behaviour is also 
reported in Ref. [8]. Another interesting feature of the curves shown in Fig. [T]is that 
for T < T c 57 MeV there is a significant decrease of the transition density, which is 
associated with the appearance of color superconductivity below T c . This results are 
in agreement with similar calculations that employ the Nambu-Jona-Lasinio model for 
quark matter [15] . As expected, we found that the transition density increases when 
the drop radius decreases, as a consequence of the larger contribution of the surface and 
curvature terms in Eq. f l33l . 

The behaviour of the transition density with the gap parameter Ao is shown in 
Fig. [2] for different values of /z^, T and B. Notice that the transition density is a 
decreasing function of the gap parameter A . The effect is strong, e.g. the transition 
density for Ao ~ 100 MeV is much smaller than for Ao = MeV. For sufficiently small 
A the transition density has constant values. This is because this part of the curve 
corresponds to temperatures that are larger than the critical temperature T c = 0.57Ao, 
and therefore the pairing gap A(T) is zero. 

As can be seen in Fig. El we obtain a < 18 MeV fm~ 2 and £ ~ 4 — 12 MeV fm -1 
which are ordinarily an order of magnitude smaller than the values obtained within the 
Nambu-Jona-Lasinio model [13]. Our values for a are also smaller than the constant 
value a = 30 MeV fm -2 assumed in Ref. pQ. 

In Figs. H] and [5] we present the results for the nucleation rate given by Eq. fl39|) . 
From Fig. H]we see that there are large variations in the value of log 10 (r[cm~ 3 s -1 ]) for 
different values of B, but the qualitative behaviour is nearly the same in all cases: (i) 
the nucleation rate increases hugely when the temperature falls from 61 MeV to 1 MeV, 
and (ii) the nucleation of droplets larger than some value R* are strongly suppressed 
(typically R* ~ 100 — 200 fm). On the contrary, droplets with R < R* have a huge 
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Figure 1. Mass-energy density and temperature of hadronic matter at which 
deconfinement is energetically favoured for R = oo and R = 2 fm. The deconfmement 
density for finite size drops is larger than for the bulk case due to surface tension and 
curvature energy. We have calculated the transition density for several radii and found 
that for R ~ 100 — 200 fm the curves are almost coincident with the case R = oo. The 
density is given in units of the nuclear saturation density po. We use Ao = 100 MeV 
and B [MeV/fm 3 ] = 100, 160 and 353. We considered different chemical potentials 
for trapped electron neutrinos in hadron matter but for simplicity we show only the 
results for /x£ [MeV] = and 200. 



nucleation rate at low enough T, and in practice should nucleate instantaneously when 
the bulk transition density is attained at the core of a protoneutron star. In Fig. [5] we 
show that the main conclusions obtained from Fig. 0] are valid for different values of 
A . Specifically we notice on the right panel that log lo r is a large negative number 
for radii above R*. On the left panel, we see that log 10 T is a large positive number for 
radii below R*, with the possible exception of some very small drops with R < 10 fm. 




Figure 2. Mass-energy density at which deconfinement is energetically favoured as 
a function of the parameter Aq. There is a strong decrease of the transition density 
p H / Po for large enough Aq. 
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Figure 3. Surface tension a and curvature energy £ for droplets with R = 2 fm and 
R = oo. The curves were calculated using the set of values of the thermodynamic 
variables that arise from Eqs. (I5T1 - I55)) . i.e. the same set of values that led to Fig. 1. 
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Figure 4. Nucleation rate of droplets as a function of their radii, for Ao = 100 MeV, 
B [MeV/fm 3 ] = 100, 160, 353, tff [MeV] = 0, 100 and T [MeV] = 1, 61. 
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Figure 5. Nucleation rate of droplets as a function of Ao for B [MeV/fm 3 ] = 160, /i^ 
[MeV] = 100 and T [MeV] = 11. 



6. Summary and Conclusions 

In the present work we analyse the nucleation of quark droplets inside a protoneutron 
star focusing our analysis in the color superconductivity, trapped neutrinos, and finite 
size effects. To describe strongly interacting matter we consider a different model for 
each phase. For the hadronic phase we use the Walecka model including the baryonic 
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octet, electrons and electron neutrinos in equilibrium under weak interactions. The just- 
deconfined quark matter is described using the MIT bag model (considering u, d and s 
flavours, electrons and electron neutrinos). To determine the energy density at which 
the deconfinement occurs we assume a first order phase transition, we impose flavour 
conservation during the transition, and we consider the deconfined phase to be in the 
energetically favoured 2SC state. The effects of finite size of the nucleated droplet are 
included through the MRE formalism, i.e. the surface tension and the curvature energy 
are not treated as free parameters but are calculated self-consistently as a function of 
the thermodynamic state of the system. 

Firstly, we determine the transition density for given values of the droplet radius, 
the gap parameter A , the temperature and the chemical potential of the neutrinos in 
the hadronic phase (see Figs. [Hand [2]). Due to the effect of color superconductivity, 
there is a significant decrease of the transition density at low temperatures. In most 
cases, our values for the surface tension (see Fig. [3]) turn out to be a factor ~ 2 — 10 
smaller that the constant value assumed in Ref. pQ within the MIT bag model and 
between one and two orders of magnitude smaller than the value calculated in [13] 
using the MRE together with the NJL model. Secondly, we calculate the nucleation 
rate of quark droplets as a function of their radii for different temperatures, pairing 
gaps A , and trapped neutrino abundances (see Figs. @]and[5]). Large droplets have 
a low transition density because surface and curvature effects tend to be small; but at 
the same time, their nucleation probability is suppressed because there is a large free 
energy involved in their formation (due to their large volume). As a consequence, the 
nucleation rate V is negligible for droplets larger than R* ~ 100 — 200 fm. We also show 
that T is huge for most radii below R* provided that the temperature is low enough. 

All these results are in qualitative agreement with previous calculations obtained 
within the frame of the Nambu-Jona-Lasinio model [13]. Quantitatively, the here 
obtained nucleation rate is very different to the one obtained within the NJL model; 
nonetheless, both are negligible above some radius R* and huge below it. Moreover, 
in spite of the large differences in a, ( and T for the MIT and the NJL models, the 
radius R* is within the same order of magnitude in both cases (some hundreds of fm). 
Since the transition density for R ~ 100 — 200 fm is almost coincident with the bulk 
transition density and the nucleation rate of such droplets is huge, the deconfinement 
transition should be triggered instantaneously when a density slightly larger than the 
bulk transition density shown in Fig. 1 is reached at some layer of a PNS. Since color 
superconductivity lowers the transition density at low temperatures, the transition is 
likely to occur after significant cooling in a massive enough protoneutron star. 
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